Aerodynamics of the Mars Microprobe Entry Vehicles 
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The selection of the unique aeroshell shape for the Mars Microprobes is discussed. A 
description of its aerodynamics in hypersonic rarefied, hypersonic continuum, supersonic 
and transonic flow regimes is then presented. This description is based on Direct Simu- 
lation Monte Carlo analyses in the rarefied-flow regime, thermochemical nonequilibrium 
Computational Fluid Dynamics in the hypersonic regime, existing wind tunnel data in the 
supersonic and transonic regime, additional computational work in the transonic regime, 
and finally, ballistic range data. The aeroshell is shown to possess the correct combina- 
tion of aerodynamic stability and drag to convert the probe’s initial tumbling attitude 
and high velocity at atmospheric- interface into the desired surface- impact orientation and 
velocity. 


Introduction 

W HEN the Maxs Surveyor mission is launched 
in February of 1999, it will transport not only 
its own lander and or biter to Mars, but two small soil 
penetrators. These two Mars Microprobes 1 axe the 
second of the Deep Space missions from the New Mil- 
lennium Program Office. Upon arriving at Maxs, the 
penetrators will be released from the Surveyor space- 
craft and begin a free fall to the surface. This paper 
focuses on the aerodynamic challenges associated with 
that Mars entry. In particular, the penetrators must 
be packaged within passive aeroshells which will carry 
their instrumented payloads safely through hypersonic 
rarefied, hypersonic continuum, supersonic, and tran- 
sonic flows to an impact with the surface at prescribed 
velocity and attitude. 

The objective of this paper is to review the decisions 
made in the selection of the aeroshells to fulfill these 
requirements, and then describe in depth the expected 
aerodynamics in each of the flight regimes. This de- 
scription is based on Direct Simulation Monte Carlo 
analyses (DSMC) in the rarefied-flow regime, ther- 
mochemical nonequilibrium Computational Fluid Dy- 
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namics (CFD) in the hypexsonic regime, existing wind 
tunnel data in the supersonic and transonic regime, 
additional CFD in the transonic regime, and finally, 
ballistic range data. 

Selection of the Aeroshell 

Selection of the aeroshell for Mars Microprobe re- 
quires consideration of the unique objectives of the 
mission. An enclosure is required which will safely de- 
liver the penetrator payload through entry to impact 
with the surface. The aeroshell must decelerate the 
vehicle during its descent to a prescribed impact ve- 
locity. It must possess sufficient stability to correct an 
initial tumbling attitude to one of forward-facing flight 
early in the trajectory, and it must maintain that ori- 
entation within small tolerances until impact. This 
attitude control must be accomplished passivly since 
there is no accomodation for active controls. Finally, 
the aeroshell must protect the payload from intense 
aerodynamic heating. The basis for this selection is 
discussed in the following paragraphs. 

The Viking 2 and Mars Pathfinder 3 landers, as well 
as the Pioneer- Venus 4 and Galileo 6 probes used spher- 
ically blunted cones for their forebody. Selection of 
the appropriate cone angle involves trade-offs between 
drag, heating, stability, and overall size. Blunter 
(i.e., larger half-angle) cones exhibit more drag for a 
given surface area as well as lower heating; sharper 
cones possess more stability. Both Pioneer- Venus and 
Galileo used 45-degree half-angle cones, but these mis- 
sions were entering atmospheres much denser than 
Mars. The two Maxs missions (Viking and Pathfinder) 
require the higher-drag properties of 70-degree sphere- 
cones to decelerate their heavy landers at sufficiently 
high altitudes to permit parachute deployment. Mars 
Microprobe is light (less than 3 kg) and must impact 
the surface with sufficient speed for penetrator inser- 
tion. In addition, Microprobe requires the highest 
possible stability. 

The impact-velocity problem can be discussed in 
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terms of the ballistic coefficient (B). 


where M is the vehicle mass, A is the frontal area, and 
Cd is the drag coefficient. Figure 1 compares the ex- 
pected impact velocity versus ballistic coefficient for 
this Mars entry from three degree-of-freedom entry 
simulations. The range of acceptable impact speeds 
is shown in the figure. For a velocity of 160 m/s, a 
ballistic coefficient near 28 kg/m 2 * is necessary. The 
current mass estimate for Microprobe is 2.4 ftp, and 
the size which can be attached to the Mars Surveyor 
spacecraft has diameter (D) of 0.325 m (A=0.08296 
m 2 ), thus the drag coefficient must be around 1.0: the 
value for a 45-degree cone. [Since this work was per- 
formed the diameter had increased to 0.35 m and the 
allowable mass has increased to 2.8 kg which still sug- 
gests a drag coefficient near 1.0.] 

The extent of spherical blunting of the nose has a 
minimal effect on the drag coefficient for a 45-degree 
half-angle cone. However, Fig. 2 reveals that in- 
creased nose bluntness (iZ n /D, where R n is the nose 
radius) decreases the maximum stagnation point heat- 
ing rate (q), since heat rate varies as the inverse square 
of the effective nose radius. Unfortunately, increased 
nose bluntness also decreases static stability, as shown 
by the decreasing negative values for the slope of the 
moment curve (C m ,a from Mach 1.65 wind tunnel mea- 
surements in Ref. 6). Selecting the appropriate extent 
of nose bluntness is a trade-off between heating and 
stability. For Microprobe, a nose radius equal to half 
of the overall vehicle’s base radius ( Rn/D = 0.25 ) is 
chosen. This is the same ratio as used in the Pioneer- 
Venus and Galileo probes. In an analogous manner, 
rounding the vehicle’s shoulders decreases the heat- 
ing at that location. However, rounding the shoulders 
decreases drag and stability. The compromise is to 
again use the Pioneer- Venus values: shoulder radius 
(R 8 ) equal to one-tenth the nose radius. It is pos- 
sible to optimize the nose and shoulder radii for the 
Microprobe mission, but these previously used ratios 
appear adequate, and their selection allows the use of 
an extensive body of existing aerodynamic data. The 
geometry of the Pioneer- Venus small probes is shown 
in Fig. 3. 

The hemispherical afterbody serves two purposes. 
First, since the vehicle’s initial attitude is tumbling, it 
may encounter the atmosphere while traveling back- 
wards. A hemispherical afterbody with its center at 
the vehicle’s center-of-gravity is not stable at this back- 
ward orientation and will foster rotation to a forward- 
facing attitude. Second, the shape has been shown to 
decrease the dynamic instability observed in blunt ve- 
hicles traversing the transonic flight regime 8 . Selection 
of the large hemispherical afterbody for Microprobe 
does not rely on previous mission’s designs, but is sim- 


ilar to the one used in the PAET experiment 7 . 

Regarding the backward stability issue, Pioneer- 
Venus, Galileo, Viking, and Mars Pathfinder were 
all hypersonically stable in either a forward or back- 
ward orientation. To avoid a backward entry, each 
were oriented nose first and spin-stabilized to assure 
a nose-first attitude at atmospheric interface. Since 
spin-stabilization is not an option for Microprobe, an 
afterbody which assures the vehicle does not trim in a 
rearward facing attitude must be chosen. 

The desired impact velocity of Mars Microprobe 
is just below Mach 1. Unfortunately, aeroshells like 
the 45-degree and 70-degree blunted cones mentioned 
above suffer a dynamic instability at small angles of 
attack in the Mach 1.0 to 2.0 range. 2 * 8 * 9 ’ 10 . That is, 
though they remain statically stable as they traverse 
that Mach range, an increase in incidence angles is 
observed. Since the instability is restricted to angles- 
of-attack less than 5 degrees, the increase in incidence 
is bounded and results in a wobbling motion. Un- 
fortunately for Microprobe, surface impact near Mach 
1.0 means this instability interferes with its require- 
ment of small angle-of-attack at impact. Both Viking 
and Mars Pathfinder deployed (or plan to deploy) 
parachutes at supersonic speeds to avoid this dynamics 
problem. Galileo unexpectedly traversed the transonic 
flight regime and flight data 10 indicate a growth in in- 
cidence angles beginning at Mach 2.0 and increasing 
to a maximum value near 15 degrees at Mach 1.0. 

The source of the dynamic instability is still de- 
bated. Sammonds 8 argued that a hemispherical af- 
terbody centered about the center-of-gravity would 
eliminate this problem since afterbody pressure forces 
would be directed through the c.g. such that asym- 
metric pressure distributions there would produce no 
pitching motions. He demonstrated this fact at Mach 
numbers up to 1.2. While it is unproven for larger 
Mach numbers, it is arguable that a hemispherical af- 
terbody will decrease the dynamic instability. 

Based on the above discussion, the forebody ge- 
ometry of the Mars Microprobe is specified to be a 
45-degree sphere-cone with R n = 0.08125 m, R 8 = 
0.008125 m, and D = 0.325m. The afterbody shape is 
a hemispherical section with radius 0.174 m which is 
centered about the expected c.g. location of the vehi- 
cle. The geometry is shown in Fig. 4. 

Trajectory 

Detailed analysis of the aerodynamics requires 
knowledge of the expected trajectory. A preliminary 
aerodynamic description of the Mars Microprobe was 
constructed using free-molecular aerodynamics, New- 
tonian aerodynamics, and Pioneer- Venus supersonic 
and transonic wind tunnel data. This description was 
then used in a preliminary six degree-of-freedom (6- 
DOF) trajectory simulation to create an estimated 
nominal trajectory. Altitude and Mach number ver- 


2 of 11 


American Institute of Aeronautics and Astronautics Paper 97-3658 



sus velocity from that trajectory are presented in Fig. 
5. Points from this trajectory in each flight regime 
were then selected for more detailed analysis. 


Table 1 Axisymmetric Results from DSMC. 


Free-Molecular Flow 

Free-molecular flow assumes there are no collisions 
between gas molecules in the flow field. The surface is 
impacted by free-stream particles which are diffusely 
reflected after full thermal accommodation. Unlike 
hypersonic-continuum aerodynamics where forces ex- 
erted on the body are primarily the integrated effect 
of surface pressures, free-molecular flow aerodynamics 
contain a significant contribution from shear stress. 

The two Mars Microprobes are to be released from 
the host spacecraft three minutes prior to atmospheric 
interface. As the tumbling spacecraft approach the 
outer reaches of the atmosphere, they will encounter 
widely spaced molecules. Surface impacts of these 
molecules will exert the first aerodynamic forces on 
the entry vehicle. Knudsen number {Kn) is the ratio 
of the quiescent gas’s mean free path to the vehicle’s 
diameter. At the outer reaches of the atmosphere, Kn 
is large. As long as it remains of order 10 or larger, 
the associated aerodynamic forces can be computed by 
free-molecular-flow methods. 

The free-molecular flow aerodynamics for the Mars 
Microprobe entry vehicles are presented in the form of 
nondimensional coefficients in Fig. 6. Since the geom- 
etry is axisymmetric, the static aerodynamics can be 
described by an axial-force coefficient (Ca), a normal- 
force coefficient (Cjv), and a moment coefficient (C m ). 
The reference area is the frontal area of the aeroshell 
(0.08296 ro 2 ), and the reference length is the diameter 
(0.325 m). Unless otherwise specified, moments are 
taken about the actual nose of the vehicle. For a c.g. 
location 0.24 D back from the nose or less, the vehicle 
is statically stable only when facing forward. 

Hypersonic Transitional and 
Continuum Aerodynamics 

As the descent continues into the upper atmosphere, 
the mean free path (and Kn) decrease and collisions 
between particles must be accounted for in describing 
forces on the vehicles. This is the transitional flow 
regime. Even lower (below 55 km altitude for Micro- 
probe), Kn drops below 0.001 and continuum methods 
can be used to describe the flow about, and the forces 
on, the vehicle. Figure 7 presents Kn versus altitude 
for the high altitude portion of the Microprobe trajec- 
tory. 

For a given Kn , the transitional aerodynamics can 
be computed using DSMC methods 11 . Axisymmetric 
solutions at eight points during the high-altitude por- 
tion of the trajectory were computed. The predictions 
for Ca are listed in Table 1. and seven of the points 
are plotted in Fig. 8. The figure also contains the free- 


t(s) 

Alt (km) 

M 

V(m/s) 

Kn 

C A 

10.4 

126.7 

33.3 

6909. 

86.5 

2.054 

20.3 

113.2 

35.1 

6916. 

10.08 

2.010 

31.8 

100.6 

37.4 

6921. 

1.25 

1.901 

38.0 

91.4 

39.1 

6923. 

0.305 

1.771 

48.0 

80.3 

37.7 

6908. 

5.46e-02 

1.536 

50.5 

73.4 

35.9 

6884. 

2.45e-02 

1.352 

65.0 

64.3 

35.6 

6751. 

6.78e-03 

1.190 

77.2 

54.9 

31.8 

6168. 

1.83e-03 

1.108 


molecular prediction as well as eight continuum CFD 
solutions discussed later. 

Values from a bridging function, selected as appro- 
priate to recreate the data when 10.0 > Kn > 0.001, 
are included in Fig. 8. The bridging function for the 
axial-coefficient is: 

Ca = C A ,c + {C A jm ~ C A ,c)sin 2 <j) ( 2 ) 


/ 3 1 

= \8 + 8 


logi Q Kn 


(7a, / m and (7a, c are the axial coefficients at the free- 
molecular and continuum limits, respectively. This 
equation approximates the monotonic decrease in Ca 
across the transitional regime. There is a region 
around Kn = 0.01 where both DSMC and continuum 
CFD methods can predict the axial coefficient. Above 
Kn of 0.01, the continuum CFD calculations become 
increasingly suspect. (Slip boundary conditions would 
be required in continuum methods to predict the flow- 
fields at such low density conditions.) 

Three-dimensional DSMC solutions over a range of 
angles-of-attack were computed at three of the points 
from Table 1. The results (for Ca , CW 5 and <7 m ) of 
these solutions are listed in Table 2 and compared with 
bridging functions similar to Eqns. 2-3 in Figs. 9-11. 
The bridging functions approximate the predictions of 
Ca , CW) and C m at angle-of-attack in this regime. 

In the hypersonic continuum regime, the Langley 
Aerothermodynamic Upwind Relaxation Algorithm 
(LAURA) CFD tool can describe the static aerody- 
namics of Microprobe. LAURA is an upwind-biased, 
point-implicit relaxation algorithm 12 for obtaining the 
numerical solution to the Navier-Stokes equations for 
three-dimensional viscous hypersonic flows in thermo- 
chemical nonequilibrium. The Mars atmosphere ver- 
sion of the code 13 contains an 8-species CO 2 — N 2 
chemical-kinetics model. This is the same computa- 
tional tool used to describe the aerodynamics of Mars 
Pathfinder 14,15 . Predictions from the code have been 
shown to be in agreement with Viking flight data 15 . 

LAURA solutions are computed at eight points in 
the trajectory’s hypersonic continuum regime. Table 
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Table 2 Three Dimensional Results from DSMC. 


Alt(km) 

a (deg) 

Ca 

C N 

C m | 

126.7 

0 

2.054 

0.000 

.0001 

126.7 

5 

2.040 

0.176 

-.0718 

126.7 

15 

1.921 

0.505 

-.2061 

126.7 

30 

1.574 

0.891 

-.3700 

126.7 

45 

1.118 

1.097 

-.4755 

126.7 

60 

0.723 

1.225 

-.5593 

126.7 

75 

0.360 

1.306 

-.5413 

126.7 

90 

.0015 

1.336 

-.5469 

126.7 

105 

-0.358 

1.305 

-.5370 

126.7 

120 

-0.725 

1.224 

-.5152 

126.7 

135 

-1.122 

1.095 

-.4672 

126.7 

150 

-1.570 

0.890 

-.3626 

126.7 

165 

-1.921 

0.506 

-.2040 

126.7 

180 

-2.054 

0.000 

-.0002 

100.6 

0 

1.901 

0.000 

-.0004 

100.6 

5 

1.891 

0.156 

-.0657 

100.6 

15 

1.788 

0.452 

-.1910 

100.6 

30 

1.470 

0.794 

-.3403 

100.6 

45 

1.071 

0.991 

-.4391 

100.6 

60 

0.716 

1.126 

-.4911 

100.6 

75 

0.369 

1.222 

-.5161 

100.6 

90 

0.008 

1.252 

-.5138 

80.3 

0 

1.536 

0.001 

-.0013 

80.3 

5 

1.520 

0.128 

-.0580 

80.3 

15 

1.424 

0.360 

-.1650 

80.3 

30 

1.180 

0.607 

-.2838 

80.3 

45 

0.904 

0.734 

-.3525 

80.3 

60 

0.656 

0.846 

-.4035 


Table 3 Axisymmetric Results from LAURA 
CFD. 


t(s) 

Alt(km) 

M 

V(m/s) 

Kn 

c A 

48 

80.5 

37.7 

6913 

6.1e-02 

1.598 

55 

73.4 

35.9 

6884 

2.5e-02 

1.378 

72 

55.8 

32.8 

6271 

2.0e-03 

1.086 

82 

48.7 

25.9 

5419 

1.0e-03 

1.061 

92 

43.3 

19.9 

4327 

6.2e-04 

1.049 

116 

34.9 

10.0 

2206 

2.6e-04 

1.027 

125 

32.6 

7.5 

1665 

2.1e-04 

1.028 

142 

29.5 

5.0 

1121 

1.6e-04 

1.035 


3 presents the zero-degree angle-of-attack Ca results. 
Additional forebody solutions are computed at 10 de- 
grees angle-of-attack for six of these trajectory points. 
Those results are listed in Table 4. 

The variation in axial coefficient at zero de- 
grees angle-of-attack across the hypersonic-continuum 
regime is presented in Fig. 12. Little change in Ca 
occurs across the Mach range. The variation at 10 
degrees angle-of-attack of the normal-force and mo- 


Table 4 3-D Results from LAURA, a — 10°. 


t(s) 

M 

Ca 

Cn 

C m 

72 

32.8 

1.0281 

0.1687 

-0.1106 

82 

25.9 

1.0353 

0.1689 

-0.1103 

92 

19.9 

1.0141 

0.1706 

-0.1125 

116 

10.0 

1.0016 

0.1691 

-0.1115 

126 

7.5 

0.9987 

0.1668 

-0.1125 

142 

5.0 

1.0133 

0.1635 

-0.1064 


ment coefficients are given in Figs. 13 and 14. Again, 
the aerodynamic coefficients vary little down to Mach 
5. (The aerodynamics below Mach 5 are discussed in 
subsequent sections.) The values predicted from the 
Newtonian flow approximation are included in Figs. 
12-14. Newtonian flow assumes that the free-stream 
flow is turned parallel to the surface. Local pressure is 
then a function only of the local surface inclination rel- 
ative to the free stream. While the simple Newtonian 
assumption agrees with the CFD precitions for the ax- 
ial and normal-force coefficients across the hypersonic 
regime, it under predicts the the absolute value of the 
moment coefficient by as much as 14 percent. 

At the higher Mach numbers, the bow shock is 
sufficiently energetic to initiate dissociation of the tri- 
atomic C0 2 molecule. The degree of CO 2 dissocia- 
tion on the stagnation streamline associated with four 
LAURA solutions is shown in Fig. 15. At Mach 33, 
C0 2 dissociates from its free-stream mass fraction of 
0.97 to 0.19 at the edge of the boundary layer. (The 
fully catalytic wall boundary condition forces complete 
recombination of the dissociated molecules.) Lesser 
degrees of dissociation occur with decreasing Mach 
number. No dissociation is predicted at Mach 10. 
When dissociation occurs, the residence time of the gas 
across the thin shock layer is insufficient to allow the 
gas to equilibrate. This phenomena is a result of the 
small length scales associated with Microprobe. There 
does not appear to be a Mach range where an equilib- 
rium gas assumption would be useful in predicting the 
flowfield about Microprobe. 

Figure 15 also shows the associated shock stand- 
off distances. The distance at Mach 33 is large due 
to rarefaction effects, the standoff distance reaches a 
minimum below Mach 26, then grows with decreasing 
speeds. 

Supersonic Aerodynamics 

In the supersonic regime, existing wind tunnel data 
can be used to describe the aerodynamics. Nichols 
and Nirengarten 6 tested 45-degree half-angle cones in 
the Jet Propulsion Laboratory (JPL) 21-inch super- 
sonic wind tunnel at Mach 3.0, 2.0, and 1.65 and in 
the JPL 20-inch hypersonic wind tunnel at Mach 5.0. 
The tests were performed in air at Reynolds number 
(based on diameter of 3.0 in. models) of 0.3 xlO 6 
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in the hypersonic tunnel and 1.0 xlO 5 6 in the super- 
sonic tunnel. The flight Reynolds numbers expected 
for Mars Microprobe in this speed range are below 0.1 
xlO 6 in the Mars CO 2 atmosphere. While the flow in 
this regime should be that of a perfect gas, the ratio 
of specific heats for air at 1.4 is higher than the Mars 
atmosphere value of 1.33. In addition, these measure- 
ments examined a model which did not include the 
hemispherical afterbody. Afterbody effects are negli- 
gible for small angles-of-attack above Mach 5 (which 
allowed the forebody-only CFD used in the hypersonic 
aerodynamics description). Below Mach 5, however, 
the afterbody’s influence increases. As shown in Figs. 
12-14, the Mach 5 measurements compare well with 
the LAURA CFD calculations. 

Figures 16-18 focus on the Mach range below 5 and 
include wind tunnel, ballistic range, and CFD data. 
The solid line in the figures labeled “database” is in- 
cluded to show trends across the Mach range. The 
figures reveal a marked change in the aerodynamics 
across the supersonic and transonic regimes. The 
change is a result of the sonic line shifting from the 
nose region to the shoulder region. At Mach 5, the 
subsonic bubble remains on the spherical nose even at 
10 degrees incidence. By Mach 2, the sonic line has 
shifted to the shoulder even at zero degrees incidence. 
At this condition, the entire forebody shock layer is 
subsonic flow. 

Pressure distributions, and aerodynamic forces are 
affected by the sonic character of the shock layer. 
When the later above the flank is supersonic, the 
pressure distributions there is flat (a characteristic 
of conical flow). When the shock layer is subsonic, 
the elliptic nature of the flow results in higher, more 
rounded, pressure distributions. The net effect on the 
resulting forces is that the axial force increases while 
normal-force and moment coefficients decrease. This 
combination results in a decrease in the static stability 
margin of the vehicle. 

If the vehicle’s c.g. does not lie on the geometric 
symmetry line, the vehicle will trim at a non-zero 
angle-of-attack. As the Mach range between 5 and 
2 is traversed, the decrease in static stability would 
introduce pitching motions for such an offset c.g. since 
the vehicle would then trim at a larger angle-of-attack. 
Since Microprobe has strict requirements on angle- 
of-attack at impact (and impact occurs shortly after 
traversing this Mach range), c.g. offset should be min- 
imized to avoid the introduction of pitching motions 
late in the trajectory. 

Transonic Aerodynamics 

Figures 16-18 also contain the transonic aerodynam- 
ics. The probe’s impact speed at just below Mach 
1, and its requirement that this impact occur while 
at small angle-of-attack demand the best possible de- 
scription of the aeroshell’s aerodynamics in this speed 


regime. Unfortunately, transonic aerodynamics are 
complex, and testing and computations are difficult 
(as well as expensive). The description of Microprobe’s 
transonic aerodynamics combines existing wind tunnel 
data, new CFD, and recent ballistic-range tests per- 
formed for Microprobe. 

Wind tunnel data 2 * 16 in the transonic regime were 
obtained in support of the Pioneer- Venus small probes 
which have the same forebody as Microprobe. Both ve- 
hicles have hemispherical afterbodies (Figs. 3 and 4), 
though Microprobe’s is larger. In Ref. 16, static aero- 
dynamics for Pioneer- Venus models were measured at 
Mach 0.2, 0.5, 0.8, 0.95, 1.1, 1.2, and 1.35 in Lang- 
ley’s 8-foot transonic wind tunnel, and Mach 1.5, 1.8, 
and 2.16 in the Langley Unitary Plan wind tunnel. 
The reference describes testing problems encountered 
at Mach numbers 1.2, 1.35, and 1.5. Figures 16-18 
exclude the data from those Mach numbers. 

Figure 16 also contains CFD results at five points 
between Mach 0.6 and Mach 2.2 using the TLNS3D 
(Thin Layer Navier Stokes 3-Dimensional) code 17 . 
These calculations axe about the Microprobe geom- 
etry for conditions coincident with the Brooks (Ref. 
16) wind tunnel data. The results presented were pre- 
ceded by an extensive grid-resolution study to resolve 
the wake region. The CFD calculations span the range 
where the Brooks study encountered problems with re- 
flected wall-shock interactions. In an effort to assess 
the usefulness of the Pioneer- Venus wind-tunnel mea- 
surements in the transonic regime, additional TLNS3D 
CFD was computed at the same five conditions for 
the Pioneer- Venus geometry. The axial-force coeffi- 
cients for the two geometries are compared in Fig. 19. 
Higher base pressures on the large Microprobe after- 
body result in a decrease (as large as 6 percent) in the 
predicted drag coefficient for the Microprobe geometry 
relative to the Pioneer- Venus geometry. 

Ballistic range tests of the Microprobe geometry 
were performed at Wright Patterson Laboratory’s Aer- 
oballistic Research Facility (ARF) at Eglin Air Force 
Base. The test were concentrated around Mach 0.8 
and 1.4. The primary goal of these tests were to estab- 
lish the dynamic stability of the vehicle (discussed in 
the next section). Static aerodynamics, however, were 
also measured and are included in Figs. 16-18. With 
the exception of (7/v, there is good agreement among 
the ballistic range measurements, the CFD, and the 
wind tunnel data. The ballistic range is accurate at 
measuring axial forces, it has difficulty resolving the 
small normal forces associated with a blunt body such 
as Microprobe. 

The Reynolds numbers based on diameter for this 
wind tunnel and ballistic range data is between 0.5 
and 3.0 xlO 6 . The expected flight conditions have 
Reynolds numbers between 50,000 to 100,000 for this 
speed range. (Additional transonic aerodynamics have 
recently been measured in the TsNIIMash U-21 wind 
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tunnel at the correct Reynolds number and will appear 
in a subsequent publication.) 

Dynamic Derivatives 

Detailed measurement of the dynamic derivatives 
for Microprobe over its entire trajectory to a level of 
accuracy necessary for precise 6-DOF simulation is dif- 
ficult given the time and cost restraints on the project. 
Microprobe must impact the surface while flying at 
small angle-of-attack. This impact will occur around 
Mach 0.8. The primary concern for dynamic stability, 
therfore, is in the transonic flight regime. Through 
use of existing data measured for the Mars Viking 
and Pioneer- Venus missions in conjunction with the 
transonic dynamic data measured in the Eglin ARF 
ballistic range, an assessment of the transonic dynamic 
stability of Microprobe can be made. 

In preparation for the Mars Viking mission, Uselton 
et al. 9 performed extensive wind tunnel measurements 
of the dynamic stability of 60- and 70-degree half- 
angle cones using a forced-oscillation technique. That 
reference concludes there is no appreciable difference 
between the dynamic response of the two cones when 
the c.g. location of the 70 degree cone is more forward 
than that of the 60 degree cone. The reference also 
observes a decrease in dynamic stability for a given 
cone angle when the c.g. is moved further from the 
nose. Based on these two observations, it is arguable 
that the 60 degree cone is a more dynamically stable 
shape than the 70-degree cone, but this increase in 
stability can be offset by an aft shift in the c.g. If this 
reasoning is extrapolated to 45-degree cones, a 45- de- 
gree cone should exhibit more dynamic stability than a 
60-degree cone. Therefore, the Uselton 60-degree cone 
data should represent a conservative description of the 
dynamic performance of Microprobe. Figures 20-22 
present the angle-of-attack variation of the dynamic 
derivatives at three Mach numbers in the transonic 
regime. Since the geometry is symmetric, the dynamic 
derivatives can be specified by a single parameter, 
{Cm, q + 

Dynamic tests in support of Pioneer- Venus were per- 
formed in the Langley 8 foot transonic wind tunnel 
using a similar forced oscillation technique 4 . The pio- 
neer Venus small probe geometry (Fig. 3) has the same 
forebody as Microprobe. The c.g. location is 0.22 D 
back from the nose for the Pioneer- Venus tests. Those 
tests were restricted to Mach numbers below 1.3 so 
representative values are included only in Figs. 20 and 
21 . 

Figures 20 and 22 also contain the Eglin ARF bal- 
listic range data for the two Mach numbers (0.8 and 
1.4) examined . The ballistic range did not distinguish 
an angle-of-attack variation in the dynamic damp- 
ing, those results (two separate measurements for each 
Mach number) axe presented as the straight dotted 
lines in the figures. 


At Mach 0.8 (Fig. 20), ail three data sources pre- 
dict dynamic stability (a negative quantity) across 
the angle-of-attack range. The figure shows that a 
45-degree cone possess a higher degree of dynamic sta- 
bility than the 60-degree cone. The Eglin ARF data 
is in good agreement with the wind tunnel data. 

The instability at small angles-of-attack for Mach 
1.2 and 1.4 is indicated by the positive values for 
(Cm f q + Cm.oc) in Figs. 21 and 22. The ARF Ballistic 
range data at Mach 1.3-1.45 did not detect a transonic 
dynamic instability for the Microprobe geometry. This 
may be due to Microprobe’s hemispherical afterbody 
which, as reported by Sammonds, can eliminate or re- 
duce the transonic dynamic instability. 

Accurate 6-DOF simulation over the entire en- 
try trajectory requires knowledge of the dynamics at 
higher speeds. The Uselton measurements are re- 
stricted to Mach numbers at 3.0 and below. For larger 
Mach numbers, other ballistic range tests on blunt 
bodies 18 do not discern significant changes with Mach 
number. Thus, Mach 3.0 values for a 60-degree cone 
could be used as representative of the higher speed per- 
formance. Past experience with 6-DOF simulations 
have also indicated that large variations in the dy- 
namic derivatives at the higher Mach numbers have 
little effect on the hypersonic-flight dynamics of the 
vehicle. 

As with the supersonic and transonic data, there 
is a large difference between the Reynolds numbers 
for the dynamic data discussed above and the ex- 
pected flight values for Microprobe. Uselton 9 exam- 
ined Reynolds number effects on dynamic stability for 
70-degree cones. No clear trend was revealed. (Ad- 
ditional dynamic windtunnel test in the TsNIIMash 
U-21 at the correct Reynolds number will be discussed 
in a subsequent publication.) 

Conclusions 

Selection of the aeroshell for Mars Microprobe in- 
volved consideration of the unique objectives of the 
mission. An enclosure was required which possessed 
the appropriate drag and stability characteristics to 
safely deliver the penetrator payload through entry 
to impact with the surface. A 45-degree half-angle 
cone with rounded nose and shoulder was shown to 
be a suitable candidate for the forebody shape. A 
hemispherical afterbody centered about the vehicle’s 
center-of-gravity location was chosen for the afterbody. 

The aerodynamics of this aeroshell are discussed in 
each of the expected flight regimes during its Mars 
entry. This description makes use of Direct Simulation 
Monte Carlo (DSMC) analyses, Computational Fluid 
Dynamics (CFD), existing wind tunnel data, and new 
ballistic range data. 

Free-molecular flow aerodynamics indicate that the 
vehicle is statically unstable backwards. If the tum- 
bling vehicle has a rearward orientation when it en- 
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counters the atmosphere, it will reorient itself to a 
forwaxd-facing attitude. 

In the high-altitude transitional flow regime, a 
bridging function is shown to approximate the vari- 
ation of the aerodynamic coefficients from their free- 
molecular values (at Kn = 10.0), to their hypersonic 
continuum values (Kn of 0.001). In addition, good 
agreement between CFD and DSMC predictions for 
axial coefficient is observed for 2.0 x 10 -3 j Kn j 
3.0 x 10“ 2 . 

In the hypersonic-continuum regime, the Newtonian 
assumption is shown to be a fair estimate for the axial 
and normal force coefficients for angles-of-attack up to 
10 degrees, but it under predicts the magnitude of the 
moment coefficient by as much as 14 percent. 

In the hypersonic regime, the small length scales as- 
sociated with the forebody shock layer are insufficient 
to allow the dissociation of CO 2 to reach completion. 
Accurate description of the species profiles requires in- 
clusion of chemical nonequilibrium effects for all cases 
in which chemical reactions occur. At Mach 10 and 
below, no CO 2 dissociation is predicted. 

Around Mach 5, the sonic line shifts from the nose to 
the shoulder resulting in an increase in axial coefficient 
and a decrease in static stability of the vehicle. 

Ballistic range tests reveal that the combination 
of 45-degree cone, forward c.g. location, and hemi- 
spherical afterbody has reduced the dynamic instabil- 
ity often suffered by large-angle cones traversing the 
transonic-flow regime. 
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Fig. 1 BalHstic coefficient required for different 
impact speeds. 



njo 


Fig. 2 Increased nose radius trade between max- 
imum stagnation-point heating rate and transonic 
static stability. 


Ptoneer-Venus Small Probe Geometiy 



Fig. 3 Pioneer- Venus small probe geometry. 



Fig. 4 Mars Microprobe aeroshell geometry. 



Fig. 5 Estimated Mars-entry trajectory. 



Fig. 6 Free- molecular flow aerodynamics. 


8 OP 11 


American Institute of Aeronautics and Astronautics Paper 97-3658 



Fig. 7 High altitude flight regimes. 



Fig. 8 Comparison of zero angle-of-attack axial- 
force coefficients in the transitional-flow regime. 



Fig. 10 Comparison of normal-force coefficients in 
the transitional-flow regime. 



Fig. 11 Comparison of moment coefficients (about 
nose) in the transitional-flow regime. 



Fig. 9 Comparison of axial-force coefficients in the Fig. 12 Variation of axial-force coefficients in the 
transitional-flow regime. hypersonic-flow regime. 
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Fig. 13 Variation of normal- force coefficients in 
the hypersonic-flow regime. 



Fig. 14 Variation of the moment coefficient in the 
hypersonic-flow regime. 



Fig. 15 Stagnation line profiles of CO 2 at 
hypersonic-flow conditions. 



Fig. 16 Comparison of axial- force coefficient pre- 
dictions in the transonic and subsonic regimes. 



Fig. IT Comparison of normal- force coefficient 
predictions in the transonic and subsonic regimes. 



Fig. 18 Comparison of the moment coefficient pre- 
dictions in the transonic and subsonic regimes. 
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Fig. 19 Comparison of the TLNS3D CFD predic- 
tions for axial-force coefficient on different after- 
bodies. 
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Fig. 20 Comparison of dynamic derivatives in the 
transonic regime (M— 0.8). 



Fig. 21 Comparison of dynamic derivatives in the 
transonic regime (M— 1.2). 
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Fig. 22 Comparison of dynamic derivatives in the 
transonic regime (M= 1.3- 1.45). 
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